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Abstract 



> 

^1 . By means of a heavy-top-quark effective Lagrangian, we calculate the three-loop 

corrections of 0{a 2 s GpMf) to the H — > 66 partial decay width of the standard- 
model Higgs boson with intermediate mass Mh <C 2Mt. We take advantage of 
a soft-Higgs theorem to construct the relevant coefficient functions. We present 
our result both in the MS and on-shell schemes of mass renormalization. The MS 
formulation turns out to be favourable with regard to the convergence behaviour. 
ft, ■ We also test a recent idea concerning the naive non-abelianization of QCD. 
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One of the central questions of elementary particle physics is whether nature makes 
use of the Higgs mechanism of spontaneous symmetry breaking to endow the particles 
with their masses. In the framework of the minimal standard model, the Higgs boson, H, 
is the missing link to be experimentally discovered in order to complete our understanding 
of mass generation. So far, all attempts to detect the Higgs boson on its mass shell have 
been in vain, with the effect that the mass range Mh < 65.2 GeV has been ruled out at 
the 95% confidence level (CL) [|TJ. However, experimental precision tests of the standard 
electroweak theory are sensitive to the Higgs boson via quantum corrections. A recent 
global fit |H has yielded M H = 149 + g 4 2 8 GeV together with a 95% CL upper bound of 
550 GeV. 

A Higgs boson with Mu <^ 135 GeV decays dominantly to bb pairs. This decay 
mode will be crucial for Higgs-boson searches at LEP 2, the Fermilab Tevatron after the 
installation of the Main Injector, a next-generation e + e~ linear collider, and a future fi + fi~ 
collider. 

The present knowledge of quantum corrections to the H — > bb partial decay width 
has recently been reviewed in Ref. [Q. At one loop, the electroweak || and QCD |J 
corrections are known for arbitrary masses. In the limit Mh <C 2M 4 , in which we are 
interested here, the terms of 0(X t ), where X t = (GpM 2 /87v 2 ^/2j , tend to be dominant. 
They arise in part from the renormalizations of the Higgs wave function and vacuum 
expectation value, which are independent of the quark flavour || . In the case of bottom, 
there is an additional non-universal 0(X t ) contribution |3J, which partly cancels the 
flavour-independent one. At two loops, the universal and bottom-specific || (D(a s X t ) 
terms are available. Furthermore, the first [[J and second [ TD | terms of the expansion in 
M 2 /M H of the 0(a 2 ) five-flavour QCD correction have been found. As for the top-quark- 
induced correction in Oi^a 2 ), the full M t dependence of the non-singlet (double-bubble) 
contribution [[□]] as well as the first four terms of the M H /M 2 expansion of the singlet 
(double-triangle) contribution fl2| have been computed. At three loops, the 0(o? s ) non- 
singlet correction is known in the massless approximation |[L3|| . Furthermore, the H — > gg 
decay width receives a 0(a s ) correction due to the bbg final state which may also be 
interpreted as a 0(a^,M H /M^) correction to the H —> bb decay width |TJ]. 

We have taken the next step by evaluating the three- loop 0(a 2 X t ) corrections to the 
H — > bb decay width. In this letter, we report the key results. A detailed description of 
our analysis together with a discussion of the 0(a 2 X t ) corrections to the less important 
H — > qq decay widths, where q = u,d, s, c, will be presented elsewhere |15| . 

We now outline our procedure. We construct a rif = 5 effective Yukawa Lagrangian, 
jCy-i by integrating out the top quark. This Lagrangian is a linear combination of 
dimension-four operators acting in QCD with rif = 5 quark flavours, while all M t de- 
pendence is contained in the coefficient functions. We then renormalize this Lagrangian 
and, by exploiting the renormalization-group (RG) invariance of the energy-momentum 
tensor, rearrange it in such a way that the renormalized operators and coefficient func- 
tions are separately independent of the renormalization scale, /i. The final result for £|? 
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exhibits the following structure: 
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where q is runs over u, d, s, c, b and the primed objects refer to QCD with n/ = 5. Here, 5 U 
is the universal correction resulting from the renormalizations of the Higgs wave function 
and vacuum expectation value. The square brackets denote the renormalized counterparts 
of the bare rif = 5 operators 
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where G 

labels bare fields and parameters. 



afiu is the colour field strength, is the covariant derivative, and the superscript 

C 2q , and C 3 



coefficient functions. Notice that 



O' 



3? 



'3q 



are the respective renormalized 
vanishes by the fermionic equation of motion, 
so that C 3g will not appear in our final result. Nevertheless, the inclusion of 0' 3q is 
indispensable for the determination of C\ and C 2q - The RG-improved formulation thus 
obtained provides a natural separation of the rif = 5 QCD corrections at scale \x = Mh 
and the top-quark- induced n/ = 6 corrections at scale // = M t , in the sense that the final 
result for the H — >• bb decay width will not contain logarithms of the type ln(M t 2 /M^) if 
the rif = 5 and rif = 6 corrections are expanded in a^(Mn) and a^(M t ), respectively. 

In contrast to the two-loop 0(a s X t ) case |§, where it was sufficient to consider just one 
term in Cy , we need to take into account three types of operators and to allow for them 
to mix under renormalization. The mixing terms are related to the 0{a 2 s ) double-triangle 
contribution considered in Ref. |12|] and extend the latter to 0(a 2 s X t ). 

Similarly to Ref. ||, we can take advantage of the higher-order formulation [16| of a 



well-known soft-Higgs theorem [I7| to simplify the calculation of the coefficient functions. 
This allows us to relate a huge number of three-loop three-point diagrams to a manageable 
number of three-loop two-point diagrams. Specifically, we have to compute 24 irreducible 
three-loop two-point diagrams for q ^ b and, in addition, 54 ones for q = b. Typical 
examples are depicted in Fig. |I[ Such a theorem is not available for the gauge interactions, 
which might explain why three-loop 0(a 2 s X t ) corrections have not yet been calculated for 
the Z — »• qq decay widths, including the important case of Z — > bb, which has recently 
attracted much attention in connection with the so-called Rb anomaly. 

From Eq. (0) we may derive an expression for the H —>■ bb decay width, appropriate 
for Mh 2M t , which accommodates all presently known corrections, including the new 
ones of 0(a 2 s X t ). It reads 



r ( H -»• bb 
Here, 
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Figure 1: Typical diagrams generating universal and non-universal 0(a 2 X t ) correc- 
tions to T [h — > bb) . Bold-faced (dashed) lines represent top quarks (Higgs or Goldstone 
bosons) . 



is the Born result including the full mass dependence. As is well known || , we may avoid 
the appearance of large logarithms of the type \n(M H /vn^) m the QCD correction, A^ CD , 
by taking m& in Eq. (|j) to be the MS mass evaluated with rif — 5 quark flavours at 
scale /i = Mh, ra^ (Mh-). Consequently, we may put m& = in A^ CD . We may proceed 
similarly with the QED correction, A^ ED , which then takes the form 

QED= 17«(M^ 

where a(u) is the MS fine-structure constant. In turn, (Mh) is then also shifted by a 
QED correction |l8j from the pole mass, Mb. A™ eak denotes the weak correction with 

Xt=0 

the leading 0(X t ) term stripped off. If we put m h = and consider the limit M H -C 2M W , 
simplifies to M 



A weak 
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where C 2 W = I - s 2 w = M^/M 2 Z and M z (M w ) is the Z-boson (W-boson) mass. A^ CD is 
the well-known QCD correction for n/ = 5 Q, 

where as = (Mh) /ft and £ is Riemann's zeta function, with values £(2) = 7r 2 /6 and 
£(3) w 1.202. This correction originates from the class of diagrams where the Higgs 
boson directly couples to the final-state bb pair. As is well known [Q, starting at 0(a 2 ), 
T(H -> bb) also receives leading contributions from the bb and bbg cuts of the double- 
triangle diagrams where the top quark circulates in one of the triangles. In the language 
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of Eq. ([]]), where the top quark only appears in the coefficient functions, this class of 
contributions is generated by the interference diagram of the operators [0[] and [0' 2b ]. 
The absorptive part of this diagram also includes a contribution from the gg cut, which 
is well known and must be subtracted. This leads to 

■? + 8C(2) 



^QCD 
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a 5 



l ln2 mg 
3 Ml 



The would-be mass singularity proportional to m 2 (m 2 /M^) cancels if the bb(g) and gg 



decay channels are combined [|12 
Equation 



is arranged in such a way that the leading M t dependence is carried by 



l + U (C 2b ) 2 -1, 



-2b, 



(9) 



where 8 U , C±, and C 2 b are defined in the context of Eq. (||). Since we use m b (Mh) in the 
rif = 5 segments of Eq. ([3]), it appears natural to also renormalize the top-quark mass 
in the nf = 6 terms according to the MS prescription. The scale-independent definition 



m 



(6) 



(lit) is singled out, since it eliminates the RG logarithms of the type ln(/i 2 /m 2 



By analogy to X t , we define x t = [GfiAI^ 712 ^)- We shall return to the M t formulation 
later on. From Ref. [OH we know that 



5 U = x t 



+ a e 
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2C(2) - 9.598 a\ 
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where a§ = a^\fit) / n . An analytic expression for 5 U , valid for // and N c arbitrary, may 
be found in Ref. |T9| . As mentioned above, C\ and C 2b are RG invariant per construction. 
Owing to the separation of the rif = 5 physics at scale \i = M H and the rif = 6 physics 
at scale fi = fi t , they are given by [[15] 

=SS( 7£)( " ,) - 7;?(M " ) 
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where /3 (n ^(/i) 



/3 



o 



(»/), 



(yu) is the Callan-Symanzik beta function, 
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(^) is the quark-mass anomalous dimension, and |15 
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Here, we have displayed the coefficient of rif = 6, for reasons which will become clear later 
on. General expressions for C% and C 2 b, valid for /i and N c arbitrary, will be provided in 
Ref. j[S). If we expand Eq. ([□]) and substitute Eq. (|12l), Eq. @ becomes 
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The appearance of the logarithms ln(/if/M#) witnesses the loss of the RG improvement 
inherent in Eq. (|TT|) . 

As a by-product, we may derive from Eq. ([!]) a formula for the H ^ gg decay width 
which includes the 0(x t ) correction. The result is 



71 



G F M%al 



[l+Xt) 



(14) 



which is in agreement with Ref. ^ffl . 

If we express Eq. flT2"D in terms of M t , we obtain 



^os 

L/ 2fe 



1 + — Ai + X t {-3 + Ae + Ai 



^ + 52C(2) + 4C(2)ln2-^C(3) 



(15) 



where A 6 = af\M t )/n, while C? s emerges from C\ by merely replacing a& and a;* with 
Aq and Xf, respectively. From Eq. (|15|) we read off that the leading 0(X t ) term receives 
the QCD correction factor (1 — 0.333 A§ — 11.219 Aq). We thus recover a pattern similar 



to the electroweak parameter Ap |2T| and the corrections 5 U , Swwh, and 5zzh HH] to 
the l + l~H, W + W~H, and ZZH vertices, respectively. In fact, the corresponding QCD 
expansions in Aq of these four observables all have negative coefficients which dramatically 
increase in magnitude as one passes from two to three loops On the other hand, if the 
top-quark mass is renormalized in the MS scheme at scale /1 = /i t , then the respective QCD 
expansions in a§ are found to have coefficients which have variant signs and nicely group 
themselves around zero |L9| . We note in passing that the study of infrared renormalons 
22fl offers a possible theoretical explanation for this observation. In the case of C2&, 
the QCD correction factor reads (1 + 2.333 a 6 + 7.032 tig). We conclude that, also in the 
case of the bbH interaction, the QCD expansion in the on-shell scheme exhibits a worse 
convergence behaviour than the one in the MS scheme. However, the difference is less 
striking than in the previous four cases. 

Finally, we would like to test Broadhurst's rule concerning the naive non-abelianization 
of QCD [p3]| . Guided by the observation that the /^-independent term of /3q emerges from 
the coefficient of rif by multiplication with —33/2, Broadhurst conjectured that this very 
relation between the n/-independent term and the coefficient of rif approximately holds 
for any observable at next-to-leading order in QCD. In Ref. []19 , this rule was applied to 



Ap, S u , 5wwh, and Szzh, and it was found that, in all four cases, the signs and orders of 
magnitude of the n/-independent terms are correctly predicted. Except for 5zzh, these 
predictions come, in fact, very close to the true values. If we multiply the coefficients of 
rif in Eqs. (|T2]) and ( |i~5l) with —33/2, we obtain —19.041 and 84.055, which has to be 
compared with the respective r^-independent terms, —28.018 and 64.223. Once again, 
the signs and orders of magnitude of the n/-independent terms, which are usually much 
harder to compute than the coefficients of rif, come out correctly. 
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